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ABSTRACT 

We present a comprehensive study of white dwarf collisions as an avenue for creating type 
la supernovae. Using a smooth particle hydrodynamics code with a 13-isotope, a-chain nuclear 
network, we examine the resulting 56 Ni yield as a function of total mass, mass ratio, and impact 
parameter. We show that several combinations of white dwarf masses and impact parameters 
are able to produce sufficient quantities of 56 Ni to be observable at cosmological distances. We 
find the 56 Ni production in double-degenerate white dwarf collisions ranges from sub-luminous to 
the super-luminous, depending on the parameters of the collision. For all mass pairs, collisions 
with small impact parameters have the highest likelihood of detonating, but 56 Ni production 
is insensitive to this parameter in high-mass combinations, which significantly increases their 
likelihood of detection. We also find that the 56 Ni dependence on total mass and mass ratio is not 
linear, with larger mass primaries producing disproportionately more 56 Ni than their lower mass 
secondary counterparts, and symmetric pairs of masses producing more 56 Ni than asymmetric 
pairs. 
Introduction 



While the preferred mechanism for type la su- 
pernovae (SNela) involves a single white dwarf 
star accreting material from a non-degenerate 
companion (Whelan & Iben 1973; Nomoto 1982; 
Hillebrandt & Niemeyer 2000), recent observa- 
tional evidence suggests a non-negligible frac- 
tion of observed SNela may derive from double- 
degenerate progenitor scenarios. Scalzo et al. (2010) 
observed the supernova SN 2007if photomet- 
rically, and assuming no host galaxy extinc- 
tion, they found 1.6±0.1M©of 56 Ni with 0.3- 
O.5M of unburned carbon and oxygen forming 
an envelope. This 56 Ni yield implies a progen- 
itor mass of 2.4±O.2M , which is well above 
the Chandrasekhar limit - the maximum mass 
for a non-rotating white dwarf (Chandrasekhar 
1931, Pfannes et al. 2010, Yoon & Langer 2004, 
Yoon & Langer 2005). It follows that two white 
dwarfs must have been involved in the event 



1 School of Earth and Space Exploration, Arizona State 
University, P.O. Box 871404, Tempe, AZ, 85287-1404 

2 Los Alamos National Laboratories, Los Alamos, NM 
87545 

3 The Joint Institute for Nuclear Astrophysics 



that produced SN 2007if, since a single white 
dwarf cannot accrete enough material to reach 
this mass without either exploding as a SNIa 
or collapsing to form a neutron star (Yoon et 
al. 2007). Furthermore, spectroscopic observa- 
tions by Tanaka et al. (2010) suggest SN 2009dc 
produced > 1.2M of 56 Ni, depending on the as- 
sumed dust absorption. This also implies a pro- 
genitor mass > 1.4M as 0.92 M of 56 Ni is the 
greatest yield a Chandrasekhar mass can produce 
(Khokhlov et al. 1993). 

Howell et al. (2006) inferred from their obser- 
vations of SN 2003fg that - 1.3M of 56 Ni was 
produced, as did Hicken et al. (2007) in their ob- 
servations of SN 2006gz. There is a growing body 
of evidence supporting double-degenerate SNela 
progenitor systems. Since any supernova arising 
from a double-degenerate progenitor scenario may 
not fit the standard templates for SNela, these 
transients must be filtered out if SNela are to re- 
main as premier cosmological tools. To that end, 
we must develop models that give clear and de- 
tectable signatures of double-degenerate SNela to 
distinguish them from standard SNela. 

Currently, models of double-degenerate progen- 
itors are split into two, dynamically different sce- 
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narios. In the first, white dwarfs in close binaries 
lose angular momenta through gravitational radi- 
ation, ultimately merging into a thermally sup- 
ported super- Chandrasekhar object or detonating 
outright (Iben & Tutukov 1984; Webbink 1984; 
Benz et al. 1989a; Yoon et al. 2007; Pakmor et 
al. 2010). In the second, two white dwarfs col- 
lide in dense stellar systems such as globular clus- 
ter cores (Raskin et al. 2009; Rosswog et al. 2009; 
Lorn-Aguilar et al. 2009; Lorn-Aguilar et al. 2010), 
where white dwarf number densities can be as high 
as w 10 4 pc -3 . This follows from conservative 
estimates for the average globular cluster mass, 
10 6 M (Brodie & Strader 2006), and for the aver- 
age globular cluster core radius, 1.5 pc (Peterson 
& King 1975), taken together with the Salpeter 
IMF (Salpeter 1955). Assuming cluster velocity 
dispersions on the order of 10 km s _1 , this allows 
for 10 — 100, z < 1 collisions per year. Obser- 
vations by Chomiuk et al. (2008) of globular clus- 
ters in the nearby SO galaxy NGC 7457 have de- 
tected what is likely to be a SNIa remnant. Given 
the difficulty in distinguishing SNela as residing 
in galaxy field stars or in globular clusters in front 
of or behind their host galaxies (Pfahl et al. 2009), 
the frequency with which these can occur warrants 
investigation. 

Numerical simulations of white dwarf collisions 
were pioneered in Benz et al. (1989b) using a 
smooth particle hydrodynamics (SPH) code. They 
concluded from their results that white dwarf col- 
lisions were of little interest as the 56 Ni yields were 
small. However, their simulations employed an ap- 
proximate equation of state for white dwarfs and 
resolutions were low, relative to what is possible 
with current computing resources. Moreover, as 
will be discussed, the infall velocities and veloc- 
ity gradients play a crucial role in the final 56 Ni 
yields. 

More recently, Raskin et al. (2009), Rosswog et 
al. (2009), and Lorn-Aguilar et al. (2010) revisited 
collisions using up-to-date SPH codes and vastly 
more particles (8xl0 5 , 2xl0 6 , and4xl0 5 , respec- 
tively). In Raskin et al. (2009), a single mass pair 
(0.6 M x2) was explored with three impact pa- 
rameters, whereas in Rosswog et al. (2009), several 
mass pairs were examined in direct, head-on colli- 
sions. Both of these papers aimed at establishing 
double- degenerate collisions as SNela progenitors, 
finding that 56 Ni is indeed produced prodigiously 



in such collisions, lending credence to their can- 
didacy. Lorn-Aguilar et al. (2010) examined one 
mass pair (0.6 M + 0.8 M ), but at a number of 
different impact parameters, ranging from those 
that resulted in direct collisions to those that re- 
sulted in eccentric binaries, aimed at establishing 
the parameters of white dwarf coalescence arising 
from collisional dynamics. 

In this paper, we revisit the three impact pa- 
rameters studied in Raskin et al. (2009) using a 
variety of mass pairings. Using 22 combinations 
of masses and impact parameters, we aim to an- 
swer five key questions; how does 56 Ni production 
depend on 

• the total mass of the system? 

• the mass ratio of the two stars? 

• the impact parameter? 

• the infall velocities of the constituent stars? 

• tidal effects? 

While the last two of these questions can be elim- 
inated with robust initial conditions, they are 
nevertheless important details that are sometimes 
overlooked. Armed with this information, we will 
be able to make some conclusions about the ob- 
servability of different combinations of collision 
parameters, and to determine whether the result- 
ing transient of any particular collision is as lumi- 
nous as a SNIa. 

The structure of this paper is as follows. In §2, 
we discuss the details of our initial conditions and 
our new hybrid burning nuclear network. In §3, 
we give the details of the results of each simulation 
that resulted in a detonation along with a study 
of the effect of numerical parameters on the 56 Ni 
yield in §3.1.2, and in §4, we discuss those that 
resulted in remnants. Finally, in §5, we summarize 
our results and conclusions. 

2. Method 

2.1. Particle Setups & Initial Conditions 

As in Raskin et al. (2009), we employ a ver- 
sion of a 3D SPH code called SNSPH (Fryer et 
al. 2006). SPH codes are particularly well suited 
to these kinds of simulations as the white dwarf 
stars involved are very dense and moving very 
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rapidly. Advecting rapidly moving, isothermal, 
cold white dwarfs in Eulerian, grid-based codes 
introduces perturbations that can be challenging 
to overcome. Moreover, because many of our sim- 
ulations are grazing impacts, conservation of an- 
gular momentum is crucial to the final outcomes, 
for which SNSPH excels (Fryer et al 2006). 

In our previous work, we used a Weighted 
Voronoi Tessellations method (WVT, Diehl & 
Statler 2006) for our particle setups. This method 
arranges particles in a pseudo-random spatial 
distribution with thermodynamic quantities that 
are consistent with the chosen equation of state 
(EOS). The default operation for this method is 
to allow the masses of particles to vary in order 
to keep their sizes, or smoothing lengths (ft), con- 
stant throughout the initial setup. This approach 
has its advantages when it comes to spatial ar- 
rangement, but one disadvantage is that it pro- 
duces a uniform level of refinement in the initial 
conditions regardless of where most of the mass 
resides. The result is that much higher particle 
counts are required to reach convergence. 

To remedy this, we modified the WVT method 
to keep mass fixed, varying h consistent with the 
density profiles of white dwarf stars. This has the 
effect of concentrating resolution where most of 
the mass resides, vastly reducing the required par- 
ticle counts for convergence. In fact, whereas in 
Raskin et al. (2009), we showed that convergence 
of the 56 Ni yield was reached at approximately 10 6 
particles, using constant mass particles, we reach 
convergence with only 200,000. A convergence test 
on particle count of our fiducial case, 0.64 M x2 
with zero impact parameter, is discussed in §3.1.2. 

A further modification we have added to our 
previous approach is an isothermalization step in 
our initial conditions. When mapping ID profiles 
for cold white dwarfs onto a resolution limited, 
3D particle setup, there is often a relaxation time, 
during which the stars oscillate before finding 
equilibrium. For white dwarf masses of « 0.6 M , 
this settling time is short, but for larger masses, 
the oscillations can continue for several minutes or 
hours. These repeated gravitational contractions 
heat the interiors of the stars until they can no 
longer be considered "cold" white dwarfs. There- 
fore, we relaxed each individual star in a modified 
version of SNSPH that artificially cools the stars 
by keeping each particle at a constant temperature 



during the stars' oscillations until they reach a cold 
equilibrium. Figure [I] shows a temperature profile 
for a 0.64 M white dwarf that has been passed 
through this isothermalization routine, indicating 
an isothermal temperature of ~ 10 T K throughout. 




x [cm] 

Fig. 1. — Temperatures and densities of particles 
lying on the x-axis in a 0.64 M white dwarf. This 
star was created using WVT and isothermalized 
to 10 7 K after ~ 5 minutes. 

As in our previous work, the initial conditions 
for the positions and velocities of the white dwarf 
stars in our simulations were generated using a 
fourth-order Runge-Kutta solver with an adaptive 
time-step that integrates simple kinematic equa- 
tions. The impactor star was initially given a small 
velocity comparable to the velocity dispersion of 
globular cluster cores, a — lOkm/s. The solver 
places the stars at O.1R apart with the proper 
velocity vectors expected for free-fall from large 
initial separations with a given velocity dispersion. 

Figure [2] compares the initial conditions of the 
0.64 M x2, head-on collision to the velocity gra- 
dient that is introduced by tidal forces. The rel- 
ative velocity of the centers of mass can be pre- 
dicted analytically for a zero impact parameter 
collision with v c = [2G(Mi +M 2 )/Ar] 1/2 , where 
Mi are the masses of the constituent white dwarfs 
and Ar is the separation of their centers of mass. 

As Figure [2] demonstrates, when the stars are 
allowed to free-fall in SNSPH from larger sep- 
arations, such as the separation of O.lR used 
throughout this paper, the velocity gradients that 
arise from tidal distortions are non-negligible. As 
will be shown in §3.1.2, the magnitudes of these 
velocities and their spreads play an important role 
in the final outcomes and the 56 Ni yields of each 
simulation as they determine how much kinetic en- 
ergy is converted to thermal energy, and thusly, 
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Fig. 2. — The velocity evolution from our initial 
conditions to the moment of first contact, indi- 
cating strong velocity gradients induced by tidal 
forces. The shaded area denotes the spread in rel- 
ative x- velocities, and v c is the relative velocity of 
the centers of mass. 



when carbon-ignition occurs. A shooting-method 
was used to determine the necessary, initial verti- 
cal separation that resulted in the final impact pa- 
rameter that we desired at the moment of impact. 
Initial velocities for all of our collision scenarios 
with zero impact parameter are given in Table [I] 
Likewise, all of our stars are initialized with 
50% 12 C and 50% 16 throughout. This approxi- 
mates typical carbon-oxygen white dwarf compo- 
sitions, and we use the Helmholtz free-energy EOS 
(Timmes & Arnett 1999; Timmes & Swesty 2000). 

2.2. Hybrid Burner 

Most large, hydrodynamic codes use some form 
of a hydrostatic nuclear network (e.g. Eggleton 
1971; Weaver et al. 1978; Arnett 1994; Fryxell et 
al. 2000; Starrfleld et al. 2000; Herwig 2004; Young 
& Arnett 2005; Nonaka et al. 2008). That is, the 
thermodynamic conditions present at the start of 
a burn calculation are not altered until the next 
hydrodynamic time step, which often times is con- 
trolled by abundance or energy changes from the 
burn calculation rather than a pure Courant con- 
dition. The effect of this is to fix the temperature- 
dependent reaction rates throughout the hydrody- 
namic time step to what they were at its start. 

There are several reasons for running a simula- 
tion this way, the most important of which is to 
avoid a decoupling of the nuclear network from the 
hydrodynamic calculation. However, for limited 
spatial, mass and time resolutions, this approxi- 
mation - that the thermodynamic conditions do 



not change rapidly enough during a burn to war- 
rant a sub-cycle recalculation of the nuclear reac- 
tion rates - fails in regimes where the nuclear re- 
actions are strongly temperature dependent, such 
as at temperatures where photo-disintegration is 
the dominant nuclear process. As Figure [3] shows 
below, the nuclear statistical equilibrium (NSE) 
state for material with p = 1 x 10 6 g cm -3 at 
Tg w 7 has most of the heavy isotopes photo- 
disintegrating back to 4 He. 
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Fig. 3. — NSE distributions for p =le7 g cm -3 
and Y e =0.5 in an a-centric nuclear network. Pro- 
ton and neutron mass fractions are plotted for ref- 
erence. At T 9 w 6, 4 He begins to dominate the 
isotope distribution. 

In such a regime, the material undergoing 
photo-disintegration experiences what amounts 
to an abrupt phase change through a strongly 
endothermic reaction. In nature, this reaction 
should rapidly cool the material before complete 
photo-disintegration, allowing these liberated de- 
part icles to react with other isotopes. However, a 
hydrostatic burn will overestimate the time-scale 
for this cooling as it assumes a full hydrodynamic 
time step is necessary for relevant pressure or 
temperature changes. With the temperature re- 
maining fixed over an artificially long time, this 
approach results in the nuclear network removing 
far too much internal energy, u, to be physical. 

Typically, one attempts to limit the impact of 
such a phase change by relying on a global time- 
step minimization scheme of the form 



At n+ i = min 



At c ,At n x f u x 



x n-l 



(i) 

where the subscript n refers to the iteration num- 
ber, At c is the Courant time, u % is the specific 
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Table 1: Initial velocities of each component star in the head-on cases of each mass pair for initial separations 
of 0.1R© . All velocities are relative to the center of mass. 



# 


mi [M ] 


m 2 [M ] 


—v\ [xl0 3 km/s] 


v 2 [xl0 3 km/s] 


1 


0.64 


0.64 


1.10 


1.10 


2 


0.64 


0.81 


1.31 


1.03 


3 


0.64 


1.06 


1.58 


0.95 


4 


0.81 


0.81 


1.24 


1.24 


5 


0.81 


1.06 


1.51 


1.15 


6 


0.96 


0.96 


1.35 


1.35 


7 


1.06 


1.06 


1.41 


1.41 


8 


0.50 


0.50 


0.97 


0.97 



internal energy of the i th particle, and f u is a di- 
mensionless parameter which constrains the maxi- 
mum allowable change in energy. Our global time- 
step is also controlled in this manner. In prac- 
tice, the conditions immediately prior to photo- 
disintegration will fix the next time-step, At n +i, 
to of order 10 _5 s for f u = 0.30. However, even 
this time-step is too large to capture the relevant 
temperature changes effecting the reaction rates 
on time-scales of order 10 _12 s. 

The alternative approach to a hydrostatic burn 
is to use a "self- heating/cooling" nuclear network 
that simultaneously integrates an energy equation 
and the abundance equation self-consistently (see 
e.g. Miiller 1986). When applied to a particle code 
like SPH, this type of calculation keeps p fixed, but 
updates temperature in a fashion that is consistent 
with the equation of state and the new internal 
energy at each sub-cycle. 

The ordinary differential equation that governs 
a hydrostatic burn calculation of the abundance 
of an isotope assuming the mass diffusion gra- 
dients are negligible, is of the form 

Yi = ^ CjRjYj 

3 

j,k J 

+ Ec^^Win*. (2) 

where the coefficients C^.j specify how many par- 
ticles of the i th species are created or destroyed, 
and are the temperature-dependent reaction 
rates for each of the different reaction types. The 
first term describes weak reactions (/3-decays and 



electron captures) and photo-disintegrations, the 
second describes two-body reactions of the type 
12 C(a,7) 16 0, and the third term describes three- 
body reactions, such as 4 He(2a,7) 12 C. The energy 
generation ODE takes the form 

e=-N A J2YimiC 2 , (3) 

i 

where mi is the rest-mass of the i th isotope (see 
e.g. Benz et al. 1989b). The density and temper- 
ature equations are simply p = and T = 0, re- 
spectively, in a hydrostatic burn. 

A self-heating at constant density calculation 
modifies only the temperature equation, starting 
from the first law of thermodynamics in specific 
mass units, 

du_P_dp = ds 

dt p 2 dt dV K ] 

where du/dt is the change in specific energy and 
ds/dt is the change in specific entropy. In ac- 
cordance with p = and employing the identity 
Ts = e, this reduces to 

dudT _ . 
dfltt 6 

f = 1 (5) 

where c v is the specific heat capacity at a constant 
volume. Equations (J2|, (J3|, and ([5| are evolved si- 
multaneously and self-consistently (Miiller 1986). 

At very high spatial resolutions and small time- 
steps, the self-heating approach would be iden- 
tical to the hydrostatic approach. As Figure [4] 
shows for the energy, temperature and composi- 
tion of a single particle over a finite and rela- 
tively large time-step as determined by Equation 
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([I]) with f u =0.30, these two burning calculations 
reach very different conclusions about the final en- 
ergy and composition of the particle after photo- 
disintegration. 




time [s] 

Fig. 4. — Calculations of the energy, temperature, 
and composition of a particle after a representa- 
tive time-step as determined by Equation ([!]) with 
f u = 0.30. Solid lines show the implicit integra- 
tions from a hydrostatic calculation, while dotted 
lines indicate those for a self-heating calculation. 
In both cases, p is kept constant, while in the self- 
heat calculation, the temperature, and thus the 
nuclear reaction rates, are recalculated at each im- 
plicit integration step, consistent with the first law 
of thermodynamics. 

To capture the relevant temperature changes 
using a hydrostatic approach would require a 
global time-step ~ 10 _12 s, where the tempera- 
tures of the two calculations have diverged by 
~ 5%. This is problematic for two reasons: 1) 
such a time-step cannot be predicted from the con- 
ditions immediately prior to photo-disintegration, 
and 2) having such a small global time-step ex- 
ceeds the limit of machine precision for many 
hydrodynamic codes, SNSPH included. When 
SNSPH attempts to calculate velocities for the 
next time-step using At ~ 10 _12 s, it often fails or 
returns zero. 

Unfortunately, a self-heating nuclear network 
can also expose the weakness of a mass resolu- 
tion limit. For a typical simulation of 10 6 parti- 
cles, each particle has a mass of ~ 10 2T g. A self- 
heating nuclear calculation for carbon-burning of 



such large masses becomes rapidly explosive on 
time-scales approaching the Courant limit. With- 
out any mechanism for energy transport on such 
short time-scales, the assumption of a homoge- 
neous burn of all 10 2T g begins to break down. The 
vigorous burning of so much material rapidly lib- 
erates more energy than the binding energy of the 
star. 

Our "middle-path" solution to these two ex- 
tremes is a hybrid-burning scheme wherein a com- 
bination of these two approaches is used under 
different circumstances. Since the hydrostatic 
approach is a better approximation for exother- 
mic reactions at our resolution limit, the self- 
heating/cooling approach is only employed for 
particles that undergo strong, net endothermic 
reactions such as photo-disintegration. This al- 
lows these particles to smoothly "step over" the 
photo-disintegration phase change without artifi- 
cially losing too much energy. We apply this ap- 
proach, along with the time-step minimization of 
Equation 0, to the a-chain aproxlS nuclear net- 
work (Timmes 1999; Timmes et al. 2000) by im- 
posing the condition that if e < after a hydro- 
static burn, the burn is recalculated employing 
Equation (J5|. 

While stepping through a photo-disintegration 
process is an interesting wrinkle for numerical sim- 
ulations of double-degenerate white dwarf colli- 
sions, it is not a significant factor for 56 Ni pro- 
duction. In all of our simulated cases, we found 
that on average < 2% of particles experienced 
this phase change. For the most part, the lo- 
cal conditions for a number of particles wherein a 
single particle might undergo photo-disintegration 
are sufficiently high-energy that the neighboring 
particles will have already initiated a detonation. 
Based on the results of our simulations, we do not 
expect that collisions with yet higher kinetic en- 
ergies than those attempted here would paradox- 
ically yield less 56 Ni due to photo-disintegration 
effects. 

3. Results & Analysis I - Detonations 

Our previous work narrowed the range of perti- 
nent impact parameters to three scenarios, which 
we revisited for each of our mass combinations. 
We simulated head-on impacts, partially grazing 
collisional impacts, and fully grazing/glancing im- 
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pacts. Table [2] summarizes the 56 Ni yields of each 
of our simulations. In this table, the impact pa- 
rameter, b (the vertical separation between the 
cores of both white dwarfs at the moment of im- 
pact) is given as the fraction of the radius of the 
primary white dwarf. Thus the 6 = column 
shows the yields for head-on impacts, the b = 1 
column indicates a full white dwarf radius and 
b = 2 indicates 2 white dwarf radii, or a fully graz- 
ing impact. 

Dursi & Timmes (2006) examined the shock- 
ignited detonation criteria for carbon in a white 
dwarf using numerical models. They derived a re- 
lationship between the density of the carbon fuel 
and the minimum radius of a burning region that 
will launch a detonation. For a carbon abundance 
of Xi2 C =0.5 and densities typically found in the 
white dwarfs used in our simulations, p ~ 10 T g 
cm -3 , their results suggest a minimum burning 
region, or match head size of ~ 10 4 cm. This is 
three orders of magnitude smaller than our small- 
est particle size, and properly resolving this crite- 
rion would require ~ 10 12 particles. Such a high- 
resolution study is too expensive with our current 
computing resources, and therefore, we acknowl- 
edge that we cannot resolve the precise detonation 
mechanism in our simulations. 

The criteria established in Dursi & Timmes 
(2006) would suggest that a single particle in any 
of our simulations can initiate a detonation. How- 
ever, in order for a detonation to be sustained, the 
energy that the initiating particle deposits in its 
neighbors must be sufficient to cause those neigh- 
bors to liberate an equal amount of nuclear energy. 
This somewhat softens the ability of a single par- 
ticle to initiate a sustainted detonation. Indeed, 
in all of our simulations, we found that several 
particles ignited nearly simultaneously, or at least 
outside of causal contact with one another in order 
to initiate a sustained detonation. Moreover, the 
pressure gradient established by particles neigh- 
boring those that reached ignition needed to be 
favorable for a significant and rapid energy depo- 
sition. 

In SPH, energy is shared between particles via 
PdV work with 

Pi 

iiij -^rrijAvij ■ ViWij, (6) 
where Pi and pi are the pressure and density of 



particle z, mj is the mass of particle j, Avij is the 
difference in velocities of particles i and j, and Wij 
is the SPH smoothing kernel. For each particle i, 
there is an implied sum over all particles j. In 
SPH formalism, the condition for a sustained det- 
onation would require that this quantity is large 
enough to ignite explosive burning in particle i. 
Put another way, if particle i generates energy 
at time to, iiji must be sufficient such that at time 
to + At, €j ~ ei, where At is the Courant time. 
This requires a proportionality between the en- 
ergy generation rate in particle i and the pressure 
gradient with its nearest neighbors, and to first 
order, this criterion reduces to 

ii > -pj^Pij. (7) 

For a given energy generation rate, large and 
positive pressure gradients can inhibit a detona- 
tion breakout. In situations where particles ig- 
nited carbon-burning, but were nevertheless un- 
able to deposit enough energy into their neigh- 
bors to cause them to also ignite, the material 
settled into a slow-burn regime rather than det- 
onating. While we cannot resolve the detonation 
mechanism to the desired precision, we compared 
one-dimensional ZND detonation profiles (see e.g. 
Fickett & Davis 1979) with detonation profiles 
from one-dimensional slices through the SNSPH 
models and concluded that our collision calcu- 
lations are resolving the detonation widths and 
structures to within 20%. 

3.1. Mass Pair 1 - 0.64 M x 2 

3.1.1. Fiducial Case 

We recalculated the « 0.6 M equal mass case 
as in Raskin et al. (2009) to establish a baseline 
comparison with our equal-mass particle config- 
uration and hybrid-burner technique. Empirical 
white dwarf mass functions {e.g. , Williams, Bolte, 
& Koester 2004) suggest collisions with this mass 
pair are expected to be among the most common. 
With our equal-mass particle constraint, the final 
mass of the star used in the simulations came to 
0.64 M . 

The right-top panel of Figure [5] shows that when 
the stars first collide the infall speeds, v x of ma- 
terial entering the shocked region are greater than 
the sounds speeds, c s , resulting in a stalled shock 
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Table 2: Simulation 56 Ni yields for various mass combinations and parameters. Values in bold are super- 
Chandrasekhar masses, and values indicated with a (D) are those simulations that resulted in remnants. 
Dashes (-) indicate combinations of parameters we did not simulate. All simulations listed here used f u = 
0.30 and 200k particles. 



# 


mi [M ] 


m 2 [M ] 


m tot [M Q ] 


b = 


b = 1 


b = 2 


1 


0.64 


0.64 


1.28 


0.51 


0.47 


U 


2 


0.64 


0.81 


1.45 


0.14 


0.53 


u 


3 


0.64 


1.06 


1.70 


0.26 


n 


u 


4 


0.81 


0.81 


1.62 


0.84 


0.84 


0.65 


5 


0.81 


1.06 


1.87 


0.90 


1.13 


n 


6 


0.96 


0.96 


1.92 


1.27 


1.32 


1.33 


7 


1.06 


1.06 


2.12 


1.71 


1.72 


1.61 


8 


0.50 


0.50 


1.00 


0.00 







in that region. The conditions in the center plane 
of this shocked region (the y-z plane, p ~ 10 6,5-7 
g cm -3 and Tg « 1) are sufficient to ignite carbon 
with an energy-generation rate scaling roughly as 
e ~ pT 22 , burning up to silicon at Tg « 3. The 
separation of material into three distinct phases is 
clearest in the left two panels of Figure |5j which 
plots particle number density in the p-T plane. 
The lower, more populated region is unshocked, 
carbon-oxygen material and is indicated in green. 
The less populated middle region, also represented 
with green at Tg « 1, is shocked material that has 
yet to reach the critical conditions for carbon ig- 
nition, and the upper, sparse region is material 
that has begun burning carbon to silicon-group 
elements, represented in red. 

The pressure gradient slopes positively in all 
directions away from the geometric center where 
this early burning begins, which is in fact at lower 
densities than the surrounding shocked medium. 
While the whole of the shocked region continues to 
heat up, causing more material to ignite near the 
center, the steep pressure gradient prevents the 
energy liberated by burning to initiate a detona- 
tion. For these burning particles, with silicon ash 
behind them and higher pressure carbon in front, 
the energy they deposit in their neighbors is in- 
sufficient to greatly alter their energy-generation 
rate. Instead, the burning region grows only as 
fast as material is heated to the critical tempera- 
tures needed for carbon ignition, Tg « 1, by the 
conversion of kinetic energy to thermal energy. 

Approximately two seconds after the stars first 
collide, sufficiently high temperatures and densi- 



ties are reached at the edges of the shocked re- 
gion to initiate carbon-burning. In these locations, 
the pressure gradient slopes negatively in all direc- 
tions. The liberated energy is free to break out, 
initiating detonations at the ignition points. The 
sound speeds in these zones are raised higher than 
the infall speeds due to the rapid rise in tempera- 
ture, and 56 Ni begins to appear in large quantities, 
indicated in blue in the left panels of Figure [6j 

Sustained detonation fronts then propagate 
through the unburned material, as well as the 
silicon "ash" that lies in the shocked region. As 
Figure [7] shows, significantly more 56 Ni is pro- 
duced during this phase. In Figure [7| it is also 
evident that the shocks overtake one another in- 
side the contact zone, shocking the material a 
second time and producing yet more 56 Ni. Less 
than one second after the detonations began, the 
entire system has become unbound, freezing out 
the nuclear reactions, as can be seen in Figure 
[8] The final 56 Ni yield for this simulation was 
0.51 M . 

In the b = 1 simulation, the added angular 
momentum distorted the shocked region between 
the two stars, resulting in detonations lighting off- 
center and off-axis as compared to the b = case. 
As Figure [9] shows, the detonation waves traveling 
through the densest portions of the shocked re- 
gions where the sound speed is highest, twist the 
material into a uniquely anisotropic configuration. 
Moreover, because much of the material is travel- 
ing nearly perpendicular to the shock, the density 
in the pre-detonation, shocked region is lower than 
in the 6 = case. This reduces 56 Ni production 
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t = 20.52770 




Density [g/cc] x [cm] 

Fig. 5. — Left Panels: Density vs. temperature for all particles from each constituent star in mass pair 
1, 0.64 M x2 and b = 0. Each point is colored indicating the isotope(s) with the greatest abundance and 
by the particle number density at each p — T coordinate. Green indicates high concentrations of carbon 
and oxygen, red indicates silicon group elements, and blue indicates iron-peak elements, most predominantly 
56 Ni. The darker the color in each group, the higher the particle number density. Right-top Panel: Sound 
speed, infall velocity, and density for particles lying on the x-axis. Right-bottom Panel: A slice in the 
x-y plane of particle densities. 
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t = 21.06990 




Density [g/cc] x [cm] 



Fig. 6. — Same format as Figure [5| at a later time in the simulation. 



t = 21.27900 




Density [g/cc] x [cm] 



Fig. 7. — Same format as Figure [5j at a later time in the simulation. 
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t = 21.68610 




Density [g/cc] x [cm] 



Fig. 8. — Same format as Figure [5j at a later time in the simulation. 



by about 7% to 0.47 M . 

Since the post-explosion, expansion phase is ho- 
mologous, the pattern of isotopes present at the 
moment the system becomes unbound is not al- 
tered by the expansion. Therefore, the velocities 
plotted in Figure [10] for several isotopes in the 
6 = and b = 1 cases are directly related to the 
radial distribution of the isotopes. 

The velocity structure preserves the isotopic 
segregation expected behind the burning front, 
with a progression from complete burning of car- 
bon & oxygen to iron-peak elements, though 
silicon-group elements, and finally ending with an 
unburned or only partially burned carbon & oxy- 
gen envelope. This layered structure is in agree- 
ment with the observations of Scalzo et al. (2010) 
and others of type la SNe suspected of having 
been produced from double-degenerate progenitor 
scenarios. 

The b = 2 scenario for this mass pair did not 
feature a detonation, and instead, resulted in a 
hot remnant embedded in a disk. Details of this 
simulation and its outcome will be discussed in §4. 



3.1.2. Variations on Parameters 

In order to assess the impact of the time-step on 
the nuclear yields, we compared three simulations 
of the 0.64 M x2, 6 = case varying the value of 
f u in Equation ([!]); one with a value of f u = 0.50, 
another with f u = 0.30, and finally, one with f u = 
0.25. As the results in table [H show for the 56 Ni 
yields, changes in the value of f u below 0.5 have 
little discernible impact on the final outcomes. 

Table 3: 56 Ni yields for 0.64 M x2, b = simu- 
lations with variations on the parameter f u and 



particle count. 


f u Particle Count 56 Ni 


0.50 


2 x 10 5 


0.51 


0.30 


2 x 10 5 


0.51 


0.25 


2 x 10 5 


0.51 


0.30 


1 x 10 4 


0.21 


0.30 


4 x 10 4 


0.31 


0.30 


4 x 10 5 


0.49 


0.30 


2 x 10 6 


0.53 



The detonations on either side of the shocked 
region are unique to the 0.64 M x2 mass pair- 



11 






Fig. 9. — A 2D slice of interpolated densities 
through the x-y plane of the 6 = 1 case of two 
0.64 M white dwarfs colliding. Four snapshots 
at different times are shown. Arrows in the top- 
left panel indicate the directions of motion of each 
star. 



ing and the 0.50 M x2 mass pairing described in 
§3.7. This is due, in large part, to the kinetic en- 
ergy at impact, which is related directly to the in- 
fall speed. With greater infall speeds, the shocked 
region heats sufficiently to initiate a detonation 
earlier, and the detonations begin nearer to the 
central region (the y-z plane). 

We tested this mechanism with a 0.64 M x2 
collision scenario by giving the constituent stars 
an artificially high infall velocity to reproduce the 
kinematic energies associated with collisions of 
larger masses. In that test, the critical temper- 
atures for carbon ignition were reached in loca- 
tions nearer to the y-z plane, but still displaced 
enough that the pressure gradient was favorable 
to a detonation. In this case, the 56 Ni production 
was actually depressed, resulting in only 0.39 M , 
due to an early detonation coupled with altered 
shock conditions. 

We also tested the effect of velocity gradients 
(tidal distortions) on the final 56 Ni yield by run- 
ning a 0.64 M x2, 6 = collision with an initial 
separation of only O.O48R with the commensu- 
rate relative velocities. In that simulation, the 
56 Ni yield was also depressed, slightly, to 0.48 M . 



i 

-1 0.1 
o 

^ 0.01 - 



0.01 



: b=o 










/ /V\\ 

/ \ V 
/ \ \v 


: b=i 






/ / A \\ 




/ ' \ \\ 1 



8.5 9 
log v [cm/s] 



9.5 



Fig. 10. — Masses of several isotopes at logarith- 
mically spaced velocity bins for the 6 = and 
6=1 cases of mass pairing 1, 0.64 M x2. 



The combination of infall velocity and tidal distor- 
tions are evidently critical for 56 Ni production. 

However, by far the most important parame- 
ter effecting the convergence of the 56 Ni yield is 
resolution. We carried out a convergence test of 
the 56 Ni yield in mass pair 1, 6 = 0, using equal- 
mass particle setups. We varied particle counts 



from 10 4 particles total, to 2 x 10 6 . As Fi gure 11 
demonstrates, convergence was reached at 2 x 10 5 
particles. This compares favorably to previous 
convergence estimates in Raskin et at. (2009) that 
concluded ~ 10 6 particles were needed for conver- 
gence using equol-h particle setups. 




Equal mass 
Equal h 



500 1000 1500 

Particle Count [xlOOO] 



2000 



Fig. 11. — Convergence of the 56 Ni yield with par- 
ticle count for simulations employing equal mass 
particles (blue, 0.64 M x 2) and equal h particles 
(red, 0.6 M x 2 from Raskin et al 2009). The 
dashed, vertical line indicates the number of par- 
ticles used in simulations throughout this paper. 
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Early work on numerical simulations of white 
dwarf collisions carried out by Benz et al (1989b) 
did not have the benefit of modern computational 
resources to reach these kinds of resolutions. Con- 
sequently, the 56 Ni yields in those simulations were 
comparatively quite low. 

3.2. Mass Pair 2 - 0.64 M + 0.81 M 

For asymmetrical collisions involving 0.64 M and 
0.81 M white dwarfs, the higher kinetic energy 
with which they collide results in several, almost 
immediate detonations near the center in the 6 = 



scenario. As Figure 12 shows, these detonation 
shocks superimpose to form a single, nearly spher- 
ical shock front that raises the sound speed above 
the infall speed for material in the 0.64 M star, 
but the pressure gradient leftward of the detona- 
tion (region 1) stalls the detonation shock, which 
can be seen as a higher-density, laminar shock at 



the rightmost edge of region 1 in Figure 12 

However, region 1 does not maintain its lentic- 
ular shape as the two stars are moving at different 
speeds relative to this shocked region. This al- 
lows the detonation shock to travel through this 
region at « Mach 1, eventually reaching fresh car- 
bon outside of region 1. This fresh carbon ignites 
explosively, creating a second detonation (region 
3 in Figure 13), which sends leading shocks back 



through region 1 and into region 2, shocking it a 
second time and eventually catching up with the 
first detonation shock. 

Most of the 56 Ni in this scenario is produced in 
the more massive star, as Table [4] demonstrates. 
Since only low-density portions of the 0.64 M star 
had entered the shocked region before the detona- 
tion, most of its contribution to the total output 
is in Si- group elements. 

In the 6 = 1 case, the pre-detonation, shocked 
region reaches much higher densities, and the 
oblique angle at which the white dwarf stars enter 
the shocked region allows more material to become 
strongly shocked by the detonation. The detona- 
tion shock also twists around the peculiar density 
contours inside the shocked region, shocking much 
of the material several times, as is seen in the 



bottom-right panel of Figure [14} The 0.81 M star 
experiences a near complete burn of all of its car- 
bon and oxygen. However, as in the 6 = case, 
most of the 0.64 M star remains unshocked at the 



Table 4: Isotope yields for the 6 = and 6 
cases of mass pairing 2, 0.64 M + 0.81 M . 



b 


Isotope 


0.81 [M ] 


0.64 [M ] 


Total [M ] 




12 Q 


0.21 


0.03 


0.24 





16q 


0.25 


0.14 


0.39 


28 Si 


0.12 


0.27 


0.39 




56 Ni 


0.13 


0.02 


0.14 




12 C 


0.02 


0.03 


0.05 


1 


16Q 


0.07 


0.14 


0.21 


28 Si 


0.12 


0.25 


0.37 




56 Ni 


0.49 


0.04 


0.53 



time of the detonation breakout. As before, the 
0.64 M star contributes mostly Si-group elements 
to the total output, as shown in Table |1J 





t=19.253799 





t=19.777100 



Fig. 14. — A 2D slice of interpolated den- 
sities through the x-y plane of the 6 = 1 
case of a 0.64 M white dwarf colliding with a 
0.81 M white dwarf. Four snapshots at different 
times are shown. Arrows in the top-left panel in- 
dicate the directions of motion of each star. 

The velocity profiles for the 6 = and 6 = 1 
cases of mass pair 2, shown in Figure [l5j rein- 
force the observation that 56 Ni is created in a 
confined region in the 6 = case, mainly in the 
densest portions of the shocked material from the 
0.81 M star. Carbon and oxygen, together, dom- 
inate the total output by mass, while in the 6=1 
case, 56 Ni is the dominant isotope, followed by 
28 Si. 
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t= 19.26010 




Density [g/cc] x [cm] 

Fig. 12. — Same format as Figure [5] for mass pair 2, 0.64 M + 0.81 M , and b = 0. Material shocked by 
the collision is labeled as region 1. Material behind the first detonation shock is labeled as region 2. 

t= 19.72830 




Density [g/cc] x [cm] 

Fig. 13. — Same format as Figure [5] for mass pair 2, 0.64 M + 0.81 M , and b = 0, at a later time in the 
simulation. Material shocked by the collision is labeled as region 1. Material behind the first detonation 
shock is labeled as region 2, and material behind the second detonation shock is labeled region 3. 
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b=0 




8.5 9 9.5 

log v [cm/s] 

Fig. 15. — Masses of several isotopes at logarith- 
mically spaced velocity bins for the 6 = and 
6 = 1 cases of mass pairing 2, 0.64 M + 0.81 M . 



3.3. Mass Pair 3 - 0.64 M + 1.06 M 

As with mass pair 2, the 6 = case of mass 
pair 3 experiences a detonation of material in the 
0.64 M star very quickly after the stars first col- 
lide. However, owing to the greater potential well 
into which the 0.64 M star is falling, the sound 
speed of the material shocked by the detonation is 
less than the infall velocity as shown in the top-left 
panel of Figure [l6j 

After « 0.7s, as the core of the 1.06 M star en- 
ters the shocked region, a second detonation lights 
on the left edge of the shocked region. This pow- 
ers a shock that travels through both stars, catch- 
ing up with the shock from the first detonation in 
the 0.64 M star. The material in the 0.64 M star 
burns mostly to 28 Si, while what burns in the 
1.06 M star burns almost entirely to 56 Ni, due to 
its higher density. The contributions from each 
star to the total elemental abundances are given 
in Table [U 



In simulations of mass pair 3 that introduced a 
non-zero impact parameter, the 1.06 M star was 
simply too compact to be significantly disrupted 
by a collision with a 0.64 M star. In both the 
6 = 1 and b = 2 cases, most of the 1.06 M star 
survived the collision, while completely disrupting 
the 0.64 M star. Details of those simulations are 
given in §4. 

3.4. Mass Pair 4 - 0.81 M x 2 

The symmetrical mass pair, 0.81 M x2 is quite 
unlike the 0.64 M x2 mass pairing discussed 
above. For the 0.81 M x2 mass pair with 6 = 0, 
several detonations occur in the y-z plane simul- 
taneously and almost immediately after impact, 
owing to the higher temperatures reached in the 
shocked region from the higher infall speeds. As 
Figure [IT] shows, these detonations superimpose 
and produce copious amounts of 56 Ni as they 
travel through the much denser material present 
inside the 0.81 M stars. This denser material 
allows for a significantly greater conversion of car- 
bon and oxygen to 56 Ni. Therefore, with only a 
26% increase in total mass of the system over the 
0.64 M x2 scenario, there is a 64% increase in 
56 Ni production to 0.84 M . 

Having denser and more compact stars also re- 
duces the sensitivity of the 56 Ni yield to impact 
parameter. Indeed, with two 0.81 M stars, both 
the 6=1 and 6 = 2 simulations resulted in detona- 
tions and significant 56 Ni production, 0.84 M and 
0.65 M , respectively. Differences in the 56 Ni 
yield for the two non-zero impact parameter sim- 
ulations stem mainly from the amount of material 
that burns to 28 Si before the detonations occur, 
with the 6 = 2 scenario featuring much more ma- 
terial burning at lower temperatures to silicon be- 
fore the detonation. The high activation energy 
of 28 Si prevents much of that material from being 
converted to 56 Ni. 



Table 5: Isotope yields for the 6 = case of mass 



pairing 3, 0.64 M + 1.06 M . 



Isotope 


1.06 [M ] 


0.64 [M ] 


Total [M ] 


12 C 


0.39 


0.02 


0.41 


16q 


0.40 


0.15 


0.55 


28 Si 


0.05 


0.24 


0.29 


56 Ni 


0.19 


0.07 


0.26 



3.5. Mass Pair 5 - 0.81 M + 1.06 M 

The 0.81 M + 1.06 M mass pair follows a very 
similar pattern to that of mass pair 2 (0.64 M + 
0.81 M ). Intermediate impact parameters allow 
more material to enter the shocked region before 
detonation, and so there is a rise in 56 Ni produc- 
tion in the 6 = 1 case over 6 = 0. However, be- 
cause both stars involved in the collision are denser 
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Fig. 16. — Same format as Figure [5] for mass pair 3, 0.64 M + 1.06 M , and b = 0. 
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Fig. 17. — Same format as Figure [5] for mass pair 4, 0.81 M x2, and 6 = 0. 



16 



than their counterparts in mass pair 2, much more 
56 Ni is produced overall. Contributions to the to- 
tal yield in the 6 = and 6 = 1 simulations are 
given m Table |6] below. 



Table 6: Isotope yields for the 6 = and b = 1 
cases of mass pairing 5, 0.81 M + 1.06 M . 



b 


Isotope 


1.06 [M ] 


0.81 [M ] 


Total [M ] 




12 C 


0.17 


0.01 


0.18 





16Q 


0.19 


0.09 


0.28 


28 Si 


0.06 


0.22 


0.28 




56 Ni 


0.58 


0.32 


0.90 




12 C 


0.05 


0.02 


0.07 


1 


16Q 


0.07 


0.09 


0.16 


28 Si 


0.06 


0.22 


0.28 




56 Ni 


0.82 


0.31 


1.13 



3.6. Mass Pairs 6 & 7 - 0. 
1.06 M x2 



,96 M x 2 & 



The O.96M x2 and 1.06 M x2 simulations 
were essentially similar to the 0.81 M x2 simula- 
tions with the exception that the greater the mass 
of the constituent stars, the less sensitive the 56 Ni 
yield was to impact parameter. Indeed, both mass 
pairs 6 and 7 produced almost the same yield in 
all three tested collision scenarios. 

What distinguishes the 1.06 M x2 mass pair 
from all the others attempted is that the 56 Ni yield 
is super-Chandrasekhar in all cases. Were such ex- 
plosions observed, there would be no doubt that 
a double-degenerate progenitor scenario of some 
kind was responsible. The resulting 1.71 M of 
56 Ni from the 1.06 M x2 simulations appear 
strikingly similar to the 1.7M of 56 Ni derived 
from the observations of Scalzo et al. (2010). 

3.7. Mass Pair 8 - 0.50 M x 2 

Finally, we studied symmetric collisions of low- 
mass, 0.50 M white dwarfs. Table [2] demonstrates 
that the 6 = collision scenario for this mass 
pairing produced less than 0.01 M of 56 Ni de- 
spite having resulted in a detonation. In this 
case, the energy generated from even mild carbon- 
burning was sufficient to unbind the stars. As in 
the 0.64 M x2, 6 = scenario, the lower veloci- 
ties with which the stars collide results in a late 
detonation. The shocked region slowly heats up 



until carbon-burning at its edges ignites a detona- 
tion. 

It is clear from the velocity profile of the most 
abundant isotopes from this collision, given in Fig- 
ure [l8j that carbon and oxygen remain mostly un- 
burned in this scenario. This seems to suggest 
that collisions of low-mass white dwarfs (M < 
0.6 M ) of the CO variety would not produce ob- 
servable transients. Other simulations introduc- 
ing impact parameters with this mass pair were 
not attempted with carbon-oxygen white dwarfs 
as the 6 = simulation yielded essentially a non- 
result. However, further investigation involving 
Helium white dwarfs is warranted. 




5.25 



7.5 8 8.5 
logv [cm/s] 



Fig. 18. — Masses of several isotopes at logarith- 
mically spaced velocity bins for the 6 = case of 
mass pairing 8, 0.50 M x2. 



4. Results & Analysis II - Remnants 

As seen in Raskin et al. (2009), the 6 = 2 
case of mass pair 1 (0.64 M x2) did not feature 
a detonation and instead formed a hot remnant 
of thermally-supported carbon and oxygen with 
some carbon-burning products. Figure [19] illus- 



trates the dynamics of this collision, starting with 
a glancing case that leads to the constituent stars 
spinning off from each other before coalescing into 
a single hot object. 

The compact remnant core after 100s featured a 
nearly constant density of p ~ 10 6 g cm -3 . It was 
surrounded by a thick, Keplerian disk w 2.0 x 10 10 
cm in radius with a scale radius ro ~ 2.3 x 10 9 
cm. The compact object at the center of the 
disk is not strictly a white dwarf since much of 
its pressure support is thermal (T w 5 x 10 8 K). 
Indeed, since degeneracy pressure support necessi- 
tates that more massive white dwarfs are smaller 
than less massive ones, this object, at « 0.8 M is 
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Fig. 19. — Snapshots of density isosurfaces at six different times for the b = 2 simulation of mass pair 1, 
0.64 M x2. After first colliding, the stars separate before coalescing into a single object. 
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far too large to be wholly degenerate (r rem « 
2.5 x 10 9 cm); larger than the 0.64 M white dwarfs 
that entered into the collision (ro.64 = 6.98 x 10 8 
cm). 

Carbon ignition nominally takes place at ap- 
proximately 7-8xl0 8 K (e.g. Gasques et al. 2007), 
but recent phenomenological models (e.g. Jiang et 
al. 2007) have suggested a strongly reduced, low- 
energy astrophysical S-factors for carbon fusion 
reactions that potentially reduce carbon ignition 
temperatures to « 3x 10 8 K, especially at densities 
of 10 9 g cm -3 . A lower carbon burning threshold 
would be of interest to future studies of collision 
remnants. 

Since the system started in a bound state, 
(T < —V, where T in this case is total kinetic 
energy and V is total gravitational potential en- 
ergy) and since any energy gained from nuclear 
processes is negligible, most of the material can- 
not escape the system and the disk remains bound 
to the compact core. It will eventually cool and 
collapse onto the surface of the compact object. 
However, the hot core may accelerate parts of the 
disk to escape velocity via radiative processes, and 
so the calculation of the final mass of the resultant 
white dwarf is beyond the scope of this paper. Suf- 
fice it to say, the final mass will not exceed the 
Chandrasekhar limit as only 1.28 M of material 
is available. 

For the b — 2 case of mass pair 2, 0.64 M + 
0.81 M , the compact remnant was slightly less 
massive at « 0.75 M . However, since the to- 
tal mass of the system is super-Chandrasekhar, 
the final remnant mass may result in a super- 
Chandrasekhar white dwarf. Again, this final 
mass will depend greatly on radiative processes, 
and the likelihood of producing a SNIa will hinge 
on the accretion rate of the disk onto the core. 

The simulations of mass pair 3, O.64M + 
1.06 , resulted in remnants in both the b = 1 
and b = 2 cases as the 1.06 M white dwarf was 
too compact for the star to be much affected by 
a grazing collision with a 0.64 M white dwarf. In 
the b = 1 case, some of the atmosphere of the 
1.06 M star was stripped away to join the mate- 
rial from the disrupted 0.64 M star in the disk, 
while in the b — 2 case, the 1.06 M star was nearly 
unaffected by the collision. Figure [20] illustrates 
that the remnant core of the b = 2 simulation 
masses « 1.0 M , and the central densities are es- 



sentially unchanged from the 1.06 M progenitor. 
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Fig. 20. — Temperature and density profiles of the 
b = 1 and b = 2 simulation remnants of mass pair 
3, 0.64 M + 1.06M©. 

Some features worth noting in the b = 2 pro- 
files in Figure [20] are the indications of a cold core 
(T core « 2 x 10 8 K) surrounded by a hot enve- 
lope (T env . % 9 x 10 8 K), and the presence of a 
strong overdensity in a part of the disk, which 
causes large spreads in density and temperature 
for the disk material. This suggests that while the 
0.64 M star is no longer a gravitationally bound 
object, it has nevertheless not been completely dis- 
rupted. 

The radial profiles of the b = 2 simulation of 
the 0.81 M + 1.06 M mass pair, shown in Figure 
[2l| indicate a very similar structure, with a high- 
temperature envelope surrounding a cold core and 
an overdensity in a part of the disk. These inho- 
mogeneities in the disk would most likely vanish 
after many crossing- times, however, it is highly 
unlikely that radiative processes or collisional ex- 
citations within the disk could remove as much 
as 0.4 M . This suggests that after a Kelvin- 
Helmholtz timescale, the remnant from the b = 2 
collision of mass pair 5 will almost certainly be 
super-Chandrasekhar. 

Some of the properties of these remnants are 
very similar to the simulated remnant explored 
in Yoon et al. (2007), which was produced by the 
merger of 0.9 M and 0.6 M white dwarfs. In that 
simulation, the remnant also featured a cold core 
(T core « 1 x 10 8 K) surrounded by a hot envelope 
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Fig. 21. — Temperature and density profiles of 
the b = 2 simulation remnant of mass pair 5, 
0.81 M +1.O6M . 

(T env . « 6 x 10 8 K), embedded inside a thick, Ke- 
plerian disk. Using a ID stellar evolution code, 
they found that such systems can indeed evolve 
on timescales ~ 10 5 yr toward becoming SNela. 

5. Discussion 

White dwarf collisions are not typically re- 
garded as SNela progenitors, and therefore, they 
have been relatively unexplored theoretically. 
Here we have conducted a comprehensive suite 
of simulations of such collisions examining the de- 
pendence of their 56 Ni yield on total mass, mass 
ratio, and impact parameter. Our results suggest 
that white dwarf collisions are a viable avenue 
for producing SNela with brightnesses that range 
from sub-luminous to super-luminous. 

In fact, in more than 75% of our simulations, 
collisions resulted in detonations, and in all but 
the least massive combination of stars, signifi- 
cant amounts of 56 Ni were produced. We found 
that even mass pairs that are below the Chan- 
draskehar limit featured explosive nuclear burn- 
ing, with the 0.64 M x2 mass pair producing 
56 Ni in quantities comparable to standard SNela. 
Moreover, the most massive combinations of stars 
produced super-luminous quantities of 56 Ni, re- 
gardless of the impact parameter, greatly increas- 
ing their likelihood of detection. The 56 Ni yields 
from these collisions are consistent with those of 
observed SNela with super- Chandrasekhar mass 
progenitors. 

Asymmetric mass pairs generally produced less 
56 Ni in head-on collisions than symmetric pairs. 
At middling impact parameters, much more 56 Ni 



was produced, however, at high impact parame- 
ters, there was little or none. This is due primarily 
to the delicate balance which must be struck be- 
tween the dynamics of the impact and the binding 
energy of the less massive star in order to establish 
a stalled shock region that can lead to a detona- 
tion. At high impact parameters, the less massive 
star is typically unbound by the collision before 
much or any 56 Ni is produced. 

For combinations of masses and impact pa- 
rameters that did not detonate, the end result 
always featured a compact, semi-degenerate ob- 
ject surrounded by a bound, thick disk of carbon 
and oxygen. Many of these systems were super- 
Chandrasekhar, and over Kelvin-Helmholtz time 
scales, these, too, are candidate progenitors for 
producing SNela. 

Our results have shown that 56 Ni production in 
white dwarf collisions is a non-linear process that 
depends on several factors, including infall veloc- 
ities and tidal distortion effects. Foremost among 
parameters to be explored in future studies is the 
composition of the constituent white dwarfs. He- 
lium has a much lower activation energy than car- 
bon or oxygen, and combinations of stars that in- 
clude helium white dwarfs would almost certainly 
produce interesting and different results. Other 
avenues to be explored include the impact of more 
detailed modeling of the isotopic profiles in the 
progenitor stars, and the possibility of sparse hy- 
drogen atmospheres. The results of these stud- 
ies will shed further light on the contribution of 
double-degenerate collisions to the observed pop- 
ulation of SNela. 
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